arXiv: 1508.05818v2 [physics.optics] 14 Oct 2016 


Excitable-like chaotic pulses in the bounded-phase regime of an opto-rf oscillator 

M. Romanelli,^’[^ A. Thorette/ M. Brunei,^ T. Erneux,^ and M. Vallet^ 

^ Institut de Physique de Rennes, Universite Rennes I - CNRS UMR 6251, Campus de Beaulieu, 35042 Rennes Cedex, Franee 
Universite Libre de Bruxelles, Optique Non Lineaire Theorique, Campus Plaine, CP 231, 1050 Bruxelles, Belgium 

(Dated: October 17, 2016) 

We report theoretical and experimental evidence of chaotic pulses with excitable-like properties 
in an opto-radiofrequency oscillator based on a self-injected dual-frequency laser. The chaotic at¬ 
tractor involved in the dynamics produces pulses that, albeit chaotic, are quite regular: They all 
have similar amplitudes, and are almost periodic in time. Thanks to these features, the system 
displays properties that are similar to those of excitable systems. In particular, the pulses exhibit a 
threshold-like response, of weh-dehned amplitude, to perturbations, and it appears possible to dehne 
a refractory time. At variance with excitability in injected lasers, here the excitable-like pulses are 
not accompanied by phase slips. 

PACS numbers: 42.65.Sf, 42.60.Mi, 42.55.-f 


I. INTRODUCTION 

Since its first appearence (in the Hodgkin-Huxley 
model of the “squid giant axon” [T]), excitability has 
proven to be a common feature of many disparate biolog¬ 
ical, chemical and physical systems, such as, for instance, 
the heart muscle, the neurons, the Belousov-Zhabotinsky 
reaction, liquid crystals, and so forth (see [IHl] and refer¬ 
ences therein). In optical systems, excitability has been 
reported for instance in optically trapped birefringent 
particles [5], optoelectronics integrated circuits [6], or in 
laser systems, in presence of optical feedback m m, of 
saturable absorbers PHD], or under optical injection El- 
lie]. The latter configuration, of particular interest for the 
present work, is conveniently described by Adler equa¬ 
tion El when the injected optical power is weak. In 
Adler model, excitability arises close to the saddle-node 
bifurcation marking the transition from phase-locking to 
phase unlocking. Therefore, the unmistakable experi¬ 
mental signature of an excitable pulse is a 27r phase 
jump dH. This is an example of an universal behavior 
because Adler equation is also well adapted to model cou¬ 
pled Josephson junctions El, biological El or microme¬ 
chanical [20] oscillators, and so forth. In the present 
work, we consider an opto-RF oscillator based on a dual¬ 
frequency laser with frequency-shifted feedback. For low 
feedback levels, the phase of this system can also be de¬ 
scribed by Adler equation, and its dynamics is analogous 
to that of optically injected lasers [21]. Another universal 
synchronization regime can be found when the coupling 
is not weak, and its effects on the amplitude of the oscil¬ 
lators cannot be neglected [22]. In this case, the transi¬ 
tion from phase-locking to phase-unlocking is not direct, 
but takes place through a window of frequency-locking 
without phase-locking, in which the relative phase is 
bounded [23]. This regime has some intriguing fea¬ 
tures |2l|25| and has attracted some interest recently. 
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not only in the context of lasers [26] [27] , but also in hy¬ 
drodynamics [28] and in nanomechanical resonators [29] . 
It is thus natural to ask whether an excitable response is 
still possible in the bounded-phase case. In the present 
work, we provide theoretical and experimental evidence 
of a mechanism leading to an excitable-like response, oc¬ 
curring at the transition from the phase-locking to the 
bounded-phase regime. Our system bifurcates from a 
phase-locked to a chaotic, self-pulsating state in which 
the pulses are not accompanied by phase slips. In “stan¬ 
dard” excitable systems, the self-pulsating state is asso¬ 
ciated to a simple attractor such as a limit cycle, and as a 
consequence the system always follows the same, unique 
path in phase space, and produces identical pulses. This 
property is not verified here. Nevertheless, albeit chaotic, 
the self-pulsating state is quite regular, because it con¬ 
sists of pulses of similar amplitudes, almost periodic in 
time (see inset of Fig.j^b) and Fig. [^(b)). Thus, we found 
that the response associated to this particular chaotic at¬ 
tractor presents features that are similar to excitability: 
Existence of a threshold i.e. need of a finite perturbation 
in order to trigger a response (which is fairly indepen¬ 
dent of the amplitude of the perturbation), and a well- 
defined refractory time during which the system cannot 
be excited again, after a first stimulus. The paper is 
organized as follows. In the next section, we describe 
the model equations and present numerical bifurcation 
curves and diagrams. We show that a chaotic, bounded 
phase attractor exists, and that close to the bifurcation 
point, intensity pulses with excitable-like characteristics 
can be found. In section [ml we describe the experimen¬ 
tal setup and results, and compare them to the numerical 
predictions. Section IV is devoted to the conclusions. 


II. MODEL AND NUMERICAL RESULTS 

A. Model equations and bifurcation diagrams 

We start our analysis from a rate-equation model that 
was introduced heuristically by Bielawski et al. [30] in 
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order to study the dynamics of a two-polarization Nd- 
doped fiber laser. The model introduces a population 
inversion for each polarization mode. Physically, this 
comes from the fact that a given active ion will interact 
preferentially with a given polarization mode (depending 
on its orientation, which is perturbed by its local envi¬ 
ronment, and, in the case of a long active medium, on 
its position along the z-direction, because of longitudi¬ 
nal spatial hole burning m)- The two populations are 
coupled by the fact that an atom in the x population 
can produce, by stimulated emission, a photon in the x 
mode, but also, with a lower probability quantified by 
a coefficient /d, a photon in the y mode. We note that 
the presence of two population inversions can be justi¬ 
fied theoretically on a more rigorous basis, starting from 
the Maxwell-Bloch equations m- This model has also 
been derived starting from a Maxwell-Bloch approach by 
Chartier et al. [32]. These rate equations permit to re¬ 
produce successfully the antiphase polarization dynam¬ 
ics of Nd- and Er-doped fiber lasers isni ES], and also 
of bulk NdiYAG lasers [34]. The antiphase oscillation 
frequency allows retrieving the value of the coupling co¬ 
efficient (3 [33|. The model equations, that wa have also 
used to analyze the synchronization dynamics of an opto- 
radiofrequency (opto-rf) oscillator, based on the beating 
between the two polarization modes of a self-injected, 
dual-frequency laser [35|, read as follows: 
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Cx^y are the amplitudes of the two laser fields coupled by 
optical injection, and the relative population inversions 
rrix^y. T] is the pump parameter, and e is the inversion life¬ 
time. The Cx field is injected in the Cy field. The injection 
process is described by two parameters, the detuning A 
and the injection strength T, which in the following are 
taken as the control parameters. The other parameters 
are P = 0.6, e = 0.0097, and y =1.2. The scaled time 
s is related to the physical time t hy s = 27r/i^t, where 
/r is the relaxation oscillation frequency. In our case 
/r 70 kHz. Physically, the Cy field is optically injected 
into Cx via an external cavity (see Fig. [^. The asso¬ 
ciated round-trip time — b ns is much smaller than 
l//i ?7 so that the coupling can be considered instanta¬ 
neous in the model. We have numerically checked the 
validity of this assumption. In Fig. we present a bi¬ 
furcation curve in the parameter plane {A, T}, computed 
using the continuation software MATCONT [36|. When 
|A| < |r|, the model admits a stable stationary so¬ 
lution, in which the laser fields are phase-locked. This 
solution bifurcates to a time-dependent state when |A| 
becomes larger than |r|. According to m, excitabil¬ 



FIG. I. Bifurcation curve. SN: saddle-node bifurcation. 
Sup (Sub) H: supercritical (subcritical) Hopf bifurcation. ZH 
and GH are codimension-two zero-Hopf and generalized-Hopf 
points respectively. 


ity arises when a system is “near a bifurcation (transi¬ 
tion) from quiescence to repetitive firing”. Therefore in 
our system excitability has to be searched close to the 
boundary of the phase-locking range. Figll] shows that 
the nature of the bifurcation to the time-dependent state 
depends on the values of A and T. For low values of A 
and T, the stationary phase-locked state is destabilized 
by a saddle-node (SN) bifurcation, leading directly to a 
phase-unlocked regime. In this range of parameters, the 
Adler approximation applies. For higher values of the 
parameters, the SN point becomes a Hopf (H) bifurca¬ 
tion point. In this case, the time-dependent state arising 
after the bifurcation corresponds to the bounded-phase 
regime, in which, contrary to the Adler case, synchroniza¬ 
tion is preserved [24] . The transition from SN to H bifur¬ 
cation occurs at the codimension-two zero-Hopf point [38] 
labelled ZH in Fig.[2 Furthermore, the Hopf bifurcation 
can be either supercritical or, for 0.36 < A ^ T < 0.97, 
subcritical. This region is of particular interest for the 
present work, because a subcritical bifurcation causes the 
system to jump to a distant region in phase space, and 
thus tends to promote a response involving large pulses. 
We underline that this subcritical bifurcation does not 
exist if the two laser modes are not coupled by the cross 
saturation P in the active medium. Indeed, for /3 = 0 
the ex field and the respective population are decoupled 
from the other variables, and the model reduces to the 
description of an optically injected laser, for which the 
bifurcation is always supercritical [35]. In this respect, 
when the feedback is not weak the system we study here 
differs fundamentally from standard optical injection. 

We now focus on the region containing the subcritical 
bifurcation, and analyze a bifurcation diagram computed 
using A as a control parameter, for a fixed value of T, that 
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FIG. 2. (a) Blue: bifurcation diagram as a function of the 
control parameter A, calculated using continuation methods. 
Red: bifurc ation diagram calculated by numerical integra¬ 
tion of eqs. ( |1|3| ). The bifurcation points are: H, Hopf; PD, 
period doubling; LP, limit point (saddle-node bifurcation); 
LPC, limit point of cycles (saddle-node bifurcation of limit cy¬ 
cles) . (b) Projection of the dynamics in a tridimensional phase 
space. The attractor is reconstructed using a time series cal¬ 
culated with a standard 4^^ order Runge-Kutta method. The 
integration step is 0.1, and the time series length is 50 000 
in normalized units. The whole time series contains about 50 
spikes. Inset: a part of the intensity time series. T = 0.9, 
and A = T + 10“^. The green and red trajectories on the 
attractor correspond to two different spikes. 


we take equal to 0.9 in the following. A rather complex 
bifurcation structure, shown in Fig. |^(a), is uncovered. 
The phase-locked state is destabilized by the subcritical 
Hopf bifurcation, where an unstable limit cycle appears. 
This cycle can be tracked by continuation, and it under¬ 
goes several secondary bifurcations without folding back 
towards the A > F region (Figj^ (a), blue curve). In¬ 
deed, numerical integration of the model equations 0i 
shows that, when A becomes larger than F, a direct tran¬ 
sition from the stationary state to deterministic chaos oc¬ 
curs [25] (Fig.|^(a), red curve). Fig.|^(b) shows a projec¬ 
tion of the chaotic attractor in the {Real(ey), Im(e^),m^} 
space, together with a time series of the beat-note inten¬ 
sity / = I p, consisting of a train of spikes that are 

quite regularly spaced. Two parts of a trajectory on the 
attractor are shown in red and green, corresponding to 
two different intensity spikes. It can be seen that, even 


if the intensity associated to the two spikes is very close, 
they actually correspond to well separated paths in phase 
space. The chaotic nature of the attractor has been as¬ 
sessed by calculating its largest Lyapunov exponent. To 
summarize, the bifurcation from quiescence to repetitive 
firing occurs at A F, from a stable steady-state (phase¬ 
locking) to a self-pulsating chaotic state. 


B. Excitable-like properties of the chaotic pulses 

The response of a system to a perturbation is called 
excitable if it presents the following features [39|. First, 
it must have a all-or-none character. This means that 
perturbations do not trigger a response if their amplitude 
is below a certain threshold, while perturbations of suf¬ 
ficient amplitude trigger a response which consists in a 
large excursion away from equilibrium, and does not de¬ 
pend on the amplitude of the perturbation. Furthermore, 
excitable systems possess a well defined refractory time 
during which they cannot be excited again, after a first 
stimulus. In standard excitable systems, the system is in 
a quiescent state close to a simple limit-cycle attractor, 
induced by a Hopf or a homoclinic bifurcation leading 
to regular oscillations or pulses. In our case, there is a 
chaotic attractor, and the pulses do not follow an unique 
path in phase space. Despite this important difference, 
the system’s response to perturbations of the detuning 
parameter A still exhibits features that are reminiscent 
of those of excitable systems (Figs.[^|^. 

First, we have checked that the response to a pertur¬ 
bation has a well-defined threshold, and is fairly indepen¬ 
dent of the perturbation’s amplitude when the threshold 
is exceeded. In order to study the response to a sin¬ 
gle perturbation, a deterministic “kick” (whose analyti¬ 
cal form is a very steep supergaussian function p{x) = 
Aexp(—x’^), with x = and n = 1000, W being the 
duration of the kick) was added to A at a given instant 
5o, and its effect on the beat-note intensity / =\^x + 
was computed (Fig. |^. Figs. (a) and Fig. Qa) show 
that the perturbation excites the relaxation oscillations. 
A large pulse is produced, followed by oscillations at the 
relaxation oscillation frequency of both the intensity and 
the population inversion. To produce the Fig. i (b), we 
have plotted the maximum of / as a function of the per¬ 
turbation amplitude A, the perturbation “energy” AW 
being held constant at the value of 1.6. For each point, 
we have integrated the equations 100 times, each time 
taking a different initial condition (the final point of the 
previous integration, i.e. a point which is very close to 
the fixed point corresponding to the locked state). This 
should mimick a real experiment repeated sequentially 
100 times. The duration of each time series was 3 000 
in normalized units. For a given value of A, some varia¬ 
tion in the response (illustrated by the error bars, visible 
only for some points of Fig. [^b)) can be observed; since 
there is no noise in the system, this variation is due to 
sensitivity to the initial conditions and is a signature of 
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the chaotic nature of the attractor. However, it is clear 
that this variation is relatively small, and that, despite 
the presence of chaos, a response of well-defined ampli¬ 
tude can be associated to a given perturbation. The am¬ 
plitude of the response shows some dependence on the 
amplitude of the perturbation, varying from around 7.8 
to 13.3 while A changes from 10“^ to 0.4. 



FIG. 3. (a) Response to a single perturbation of the detuning 
A. Upper (blue) trace: beat-note intensity. Lower (green) 
trace: population inversion, (b) Amplitude of the response as 
a function of the perturbation amplitude A. 




FIG. 4. (a) Response to two perturbations separated by a 
delay D. (b) Amplitude of the response to the second pertur¬ 
bation as a function of D. Dashed black curve: single real¬ 
ization. Red line: average over 40 realizations, ta- absolute 
refractory period, tr: relative refractory period. 

Second, we have found that, after a first excitation, 
well-defined time intervals during which the system’s re¬ 
sponse is either completely or partially inihibited can be 
identified (Fig. [^. These time intervals are reminiscent 
of the absolute and relative refractory times of excitable 
systems, as we will discuss below. In order to inves¬ 
tigate this property, we have submitted the system to 
two ’’kicks”, with variable delay D between them (see 
Fig. I^a)). Again, the crucial point is the repeatability 
of the response for a fixed delay D. So, it is mandatory 
to repeat the same sequence of two pulses many times, 
starting from different initial conditions. This is how the 
reponse curve in Fig. ib) has been obtained. In order to 
mimick a real experiment, we have taken, for each given 
value of the delay a sequence of 40 double pulses. 
Two consecutive double excitations are separated by a 
time interval of 2 000 in normalized units. The response 
as a function of averaged over 40 double excitations, 
is displayed in Figj^b). An absolute refractory period 


ta — 280 can be unambiguously identified. If D < Ta, 
then the system can never be triggered by the second 
perturbation. The absolute refractory period ta can be 
understood as follows. Looking at Fig.j^a), it can be ob¬ 
served that the intensity spike appears with a substantial 
time-delay after the kick. It is this time-delay which de¬ 
termines Ta- If the second kick arrives before that the 
intensity spike has developed, then it will not generate a 
response. The time-delay depends mainly on the distance 
from the bifurcation point after the perturbation, i.e. on 
the initial value of the detuning A and on the size of the 
perturbation, as can be observed by comparing the re¬ 
sponses in Fig.j^a) and Fig.j^a). In particular, a larger 
perturbation determines a faster response. So, also ta 
depends on the distance from the bifurcation point, i.e. 
on the detuning A and on the size of the perturbation. 
On the contrary, it does not depend much on which path 
the system follows in phase space: This makes it possi¬ 
ble to observe even in our case a feature similar to the 
refractory time of excitable systems. 

If the delay D between the two kicks is larger than 
620, a second excitation always triggers a response. In 
excitable systems, there exists also a relative refractory 
period during which the response of the system is 
weaker, but not completely inhibited m- Interestingly, 
the delay interval 280 < D < 620 can be identified, in 
a sense which will be precised in a moment, with the 
relative refractory period. The dashed black curve in 
Fig. I^b) represents the response to a single double per¬ 
turbation, as a function of D. It can be seen that, in 
the relative refractory period, the system’s response de¬ 
pends on T) in a seemingly random fashion. Indeed, we 
have found that, for a given delay inside the relative re¬ 
fractory period the response is not always the same 
for two consecutive double excitations, i.e. for different 
initial conditions. Since our model is deterministic, this 
is, again, an effect of the sensitivity of chaotic systems 
to the initial conditions. However, it can be seen that, 
when averaging over several realizations, the probability 
of triggering a second pulse increases linearly from 0 to 
1 when D goes from 280 to 620, thus producing a re¬ 
sponse curve which is very similar to the one presented, 
for instance, in Fig. 3(b) of [10], where the relative 
refractory period of an excitable semiconductor laser is 
investigated. In this sense, the [280,620] interval is rem¬ 
iniscent of the relative refractory period. We stress that 
the persistence of these features, typical of excitability, 
even in a chaotic situation, depends on the particular 
nature of the chaotic attractor we deal with. First, the 
pulses, albeit chaotic, all have very similar amplitudes, 
so that a well-defined response is obtained. Second, even 
if the pulses are chaotic, and are thus associated with 
different trajectories in phase space, a well-defined re¬ 
fractory time can still be defined, because it is essentially 
determined by the distance from the bifurcation point. 
Thanks to these characteristics, the system’s behavior is 
effectively close to an excitable one. Of course, in gen¬ 
eral this is not expected to be the case when chaos is 
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present uni- Loosely speaking, it can be expected that a 
behavior similar to excitability may be obtained only if, 
as in our case, the chaotic self-pulsating state does not 
differ too much from a perfectly periodic one. 


III. EXPERIMENTAL RESULTS 
A. Experimental setup 


(a) NdiYAG E,{y,) E/Vy) 

Pump (0.808 pm) 


Mr 


aalon ^QWP-s^M, 


Output (1.06pm) 



(Master oscillator) P-Q Analyzer 


FIG. 5. (a) Dual-Frequency Laser. Mi, 2 : cavity mirrors. 

QWP: quarter-wave plate, (b) Experimental setup that al¬ 
lows synchronizing the DFL beat-note frequency to a RF 
synthesizer. DFL: Dual-frequency laser (slave oscillator). M: 
feedback mirror. QWP: quarter-wave plate. AOM: acousto¬ 
optic modulator. P: polarizer. D: detector. P-Q analyzer: 
digital vector signal analyzer, permitting to measure the 
quadratures of I with respect to the reference signal deliv¬ 
ered by the RF synthesizer. 

The numerical results can be compared to experiments 
using the opto-radiofrequency oscillator described for in¬ 
stance in [24] . This system produces an optically-carried 
radiofrequency signal, thanks to the interference between 
the two modes of a dual-frequency laser (see Fig. [^. 
The experimental setup is as follows. The laser cavity, 
of length L = 75 mm, is closed on one side by a high- 
reflection plane mirror, coated on the 5-mm long NdiYAG 
active medium, and on the other side by a concave mirror 
(radius of curvature of 100 mm, intensity transmission of 
1% at the lasing wavelength A = 1064 nm). The active 
medium is pumped by a laser diode emitting at 808 nm. 
A 1 mm-thick silica etalon ensures single longitudinal 
mode oscillation. Two eigenmodes and Ey, polar¬ 
ized along X and with eigenfrequencies Ux and Uy re¬ 
spectively, oscillate simultaneously. An intracavity bire- 
fringent element (here two quarter-wave plates QWPs) 
induces a frequency difference, finely tunable from 0 to 
^ = 1 GHz by rotating one QWP with respect to the 
other m- Here I'y — I'x — 180 MHz. The typical output 
power of the two-frequency laser is 10 mW when pumped 
with 500 mW. When the laser output is detected by a 
photodiode after a polarizer at 45°, an electrical signal 


oscillating at the frequency difference Az^o = ~ is 

obtained. The DFL can thus be seen as an opto-RF os¬ 
cillator. In order to lock this oscillator to an external 
reference signal, we use optical frequency-shifted feed¬ 
back [42] (Fig. [^b)). The feedback cavity contains an 
acousto-optic modulator (AOM), driven by a stable RF 
synthesizer, which provides an external phase reference. 
Next, a quarter-wave plate at 45° followed by a mirror 
flips the X and y polarizations, and finally the laser beam 
is reinjected in the laser cavity after crossing again the 
AOM. As a result, a x-polarized field oscillating at the 
frequency Uy + 2/^o and a ^-polarized field oscillating at 
the frequency Vx + ‘^Jao are reinjected in the laser. We 
choose the value of 2/^o so that Vx + is close to 

Vy. Under suitable feedback conditions, Uy locks to the 
injected beam frequency z^cc+2/y^O: i-e. the frequency dif¬ 
ference Vy — Vx locks to 2/^0- We note that the optical 
reinjection has no direct effect on Ex^ because the fre¬ 
quency difference between Ux and Uy -|- 2/^o is too large. 
For the same reason, multiple round trips in the feedback 
cavity have no effect on the dynamics. The laser output 
is detected with a photodiode (3 GHz analog bandwidth) 
after a crossed polarizer, thus providing an electrical sig¬ 
nal proportional to / = \Ex Ey\‘^. The signal is then 
analyzed with an electrical spectrum analyzer, a digital 
P-Q signal analyzer and an oscilloscope. 


B. Bounded-phase pulses 

We call Adler frequency /a the maximum value of the 
detuning Au for which phase-locking can be achieved. 
In the experiments, the feedback strength is set in order 
to have /a a little smaller than the relaxation oscilla¬ 
tion frequency /r 70 kHz. Typically /a = 0.9/i? i.e. 
F = 0.9. The detuning is set in order to put the opto-rf 
oscillator very close to the boundary of the phase-locking 
range, i.e. to the Hopf bifurcation point: Au < /^. 
Typically /a — Au 100 Hz, i.e. |A — F| ^ 1. An 
experimental time series of the beat-note intensity and 
phase under these experimental condition is shown in 
Fig-E Large spiking events appear in the time series 
at random instants. When observed at a shorter time 
scale, each isolated event consists of a bunch of pulses, 
originated by strong oscillations at the relaxation oscilla¬ 
tion frequency. From the experimental phase time series, 
it is clearly seen that the phase of the beat-note signal 
is weakly affected by an intensity spike. In particular, 
the phase remains bounded throughout all the bunch of 
pulses. This differs strongly from previous reports on ex¬ 
citability generated by saddle-node bifurcations and ex¬ 
plained by Adler mechanism [T51[Til[TC]. and also from 
multipulse excitability nani], where an excitable pulse 
is necessarily accompanied by a 27r phase jump, as exper¬ 
imentally demonstrated in El- By comparing the exper¬ 
imental intensity time series and the inset of Fig. i(b), 
it appears clearly that frequency noise has to be included 
in the model in order to reproduce the experimental ob- 
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FIG. 6. Experimental intensity (top) and phase (bottom) 
time series. Note that in this figure the DC part of the signal 
has been removed for technical reasons. 




FIG. 7. Computed beat-note inte nsity (top) and phase (bot¬ 
tom) time series, using equations ( |H3| ). The parameter values 
are ^ — 0.6, fn = 70 kHz, e = 0.0097, rj — 1.2, F = 0.9, and 
A = F — 10“^ + ^(s), where ^(s) is a normally distributed 
stochastic process with a standard deviation cr = 5 10“^. 


servations. Indeed, the time interval between two chaotic 
spikes is very regular in a purely deterministic model, in 
evident contrast with the experiments. So, we interpret 
the experimental observations as follows: The oscillator 
is actually inside the phase-locking range, but very close 
to the bifurcation point, so that we observe noise-induced 
pulses as in m- A simulation taking into account these 
observations is presented in Fig. which shows a calcu¬ 
lated time series of the beat-note intensity I = 
and of the relative phase. In this simulation, the system is 
close to the boundary of the phase-locking range, and the 
detuning parameter includes an additive stochastic con¬ 
tribution ^{s) in order to account for the experimentally 


observed fluctuations of the beat-note frequency when 
the laser is free running. These fluctuations are typically 
of the order of some tens of Hz on a 1 s time scale. A 
good qualitative agreement with experiments is found. 
In particular it is confirmed that the phase of the beat- 
note signal is weakly affected by an intensity spike. The 
size of the phase excursion during the spike appears very 
close to the experimental observations. 


C. Experimental tests of excitable-like properties: 
all-or-none response, refractory period 



Time (s) 

FIG. 8. Experimental beat-note intensity as a function of the 
amplitude A of a kick to the detuning parameter. 

We have tested experimentally the excitable-like char¬ 
acter of the observed dynamics, by applying abrupt, de¬ 
terministic perturbations to the driving frequency of the 
AOM, which amounts to switching the detuning between 
two different values as in Fig. (a). For Fig. the 
system was prepared in the quiescent state, close to the 
bifurcation point, and then perturbations of increasing 
amplitude A were applied. The response is plotted as a 
function of A. When A is smaller than a certain value 
At/i, there is no response. When A is sufficiently large, 
pulses of similar amplitude are emitted, irrespective of 
the precise value of A. The existence of a threshold value 
for A, and the all-or-none character of the response ap¬ 
pear in Fig. It can also be seen that Ath varies from a 
ramp to the other. We attribute this to the fact that Afh 
is determined by the distance to the bifurcation point i.e. 
by the value of A, and in the experiment this parame¬ 
ter has inevitably some fluctuations, as explained in the 
discussion of Fig. We have also investigated the issue 
of the response to two consecutives excitations, to see if 
a property analogous to the refractory time of excitable 
systems could be observed. It was not possible to obtain 
a meaningful experimental curve to be compared to the 
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FIG. 9. Experimental response to two perturbations of the 
detuning parameter A, separated by a different time interval. 

calculated one in Fig.[^(b). Again, the reason is that the 
technical fluctuations of A have an important impact on 
the value of the absolute refractory time, expected from 
the theory to depend on the precise value of the distance 
from the bifurcation point. However, we were at least 
able to verify that, when the response to the first per¬ 
turbation happens to develop when, or just before, the 
second kick is applied, then there is systematically no re¬ 
sponse to it (Fig. (a)). On the contrary for sufficient 


delay there are two nearly identical responses. This is 
coherent with the interpretation discussed in section [HI 
and suggests that the proposed analogy with excitable 
systems is meaningful. 


IV. CONCLUSIONS 

In conclusion, we have provided experimental and nu¬ 
merical evidence of a behavior having several unusual 
features, and some properties that are reminiscent of ex¬ 
citable systems. We have observed large intensity spikes 
in the output of a driven opto-RF oscillator, and shown 
by measuring the phase variations that these events oc¬ 
cur in the bounded-phase regime. We have interpreted 
our observations as noise-induced chaotic pulses. Numer¬ 
ical calculations indicate that, despite the presence of a 
chaotic attractor, the self-pulsating state associated to it 
is sufficiently regular to produce a response with proper¬ 
ties resembling to those of excitable systems. The exper¬ 
imental results are coherent with the numerical findings, 
and suggests that the above picture is meaningful. These 
results illustrate the robustness of excitable-like proper¬ 
ties as a generic feature of nonlinear systems, since they 
can appear, as in our case, in presence of higher dimen¬ 
sionality, and even chaos. 
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